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Abstract. We give a geometric criterion for the breakdown of an Einstein 
vacuum space-time foliated by a constant mean curvature, or maximal, foli- 
ation. More precisely we show that the foliated space-time can be extended 
as long as the the second fundamental form and the first derivatives of the 
logarithm of the lapse of the foliation remain uniformly bounded. We make 
no restrictions on the size of the initial data. 



1. Introduction 

This paper is concerned with the problem of a geometric criterion for breakdown 
of solutions (M. g) of the vacuum Einstein equations. 

iW(g) - o. (i) 

To describe the problem we assume that a part of space-time M* C M is foliated 
by the level hypersurface of a time function t, monotonically increasing towards 
future, with lapse n and second fundamental form k defined by, 

k(X,Y) = - E (D X T,Y), n=(-g(Dt,Bt)y 1/2 (2) 

where T is the future unit normal to St, D is the space-time covariant derivative 
associated with g, and X, Y are tangent to St and . Let So be a fixed leaf of the t 
foliation, corresponding to t = irj. We shall refer to So as initial slice. We assume 
that the space-time region .M* is globally hyperbolic, i.e. every causal curve from a 
point p S M.* intersects So at precisely one point. We also assume that the initial 
slice verifies the following assumption. 

A 1. There exists a finite covering of So by a finite number of charts U such that 
for any fixed chart, the induced metric g verifies 

Ao^d 2 < 9ij(xMj < Aoiei 2 , VxeC/ (3) 
with Ao a fixed positive number. 

We consider the following two situations: 
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(1) The surfaces S( are asymptotically flat and maximal. 

trfc = 0. 

(2) The surfaces E t are compact, of Yamabe type —1, and of constant, negative 
mean curvature. They form what is called a (CMC) foliation . 

trk = t, t < 

Though our methods apply equally well to both situations we shall only consider 
here the latter case, which is somewhat easier to treat due to the compactness of 
the level surfaces E t . We shall thus assume in what follows that the region A4* is 
equal to Li te [ to _ t ^ t , with < 0. We can also assume that the initial hypersurfacc 
S corresponds to to = — 1. 

Remark. In the second case the CMC conjecture asserts that it should be possible 
to extend the foliation, in a smooth manner, to its maximal allowed value trfc = 
t = 0, see [And] and references therein. 

Given p 6 M* we can define a point-wise norm |II(p)| of any space-time tensor IT 
via decomposition 

X = -X°T + X, XeTM, IeTS t 
We denote by the L p norm of II on E t . More precisely, 

\\n(t)\\ LP = f \n\Pdv g 

with dv g the volume element of the metric g of £ t . The main result of this paper 
is the following theorem. 

Theorem 1.1 (Main theorem). Let (M,g) be a globally hyperbolic development of 
So foliated by the CMC level hypersurfaces of a time function t < 0, such that So 
corresponds to the level surface t = to- Assume that So verifies Al. Then the first 
time T* < 0, with respect to the t-foliation, of a breakdown is characterized by the 
condition 

lim sup (||fc(t)|| L « + ||Vlogn(i)|| L oo) = oo (4) 

More precisely the space-time together with the foliation S t can be extended beyond 
any value t* < for which, 

sup ||fc(t)|| L oo + ||VIogn(t)|| L oo =A <oo (5) 
te[t ,t*) 

Condition (5) can be reformulated in terms of the deformation tensor of the future 
unit normal T, IT = ( T )tT = £Tg- By a simple calculation, expressed relative to an 
orthonormal frame eo = T, e\, e2, e$, we find, 

7T o = 0, 7T i = n _1 Vjn, mj = -2kij. (6) 

Consistent with the statement of the main theorem we assume that T is an ap- 
proximate Killing vcctorfield in the following sense, 
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A2. There exists a constant Aq such that, 



sup ||7r(t)||£o. < A 



(7) 



te[t ,t«) 



In addition to the constant A in Al, A2 we introduce another constant IZo which 
plays an important role in the proof, which bounds the L 2 norm of the spacetime 
curvature tensor R on S , 



To prove our main theorem we have to show that if assumptions Al and A2 are 
satisfied then the space-time A4* can be extended beyond t*. We want to emphasize 
that theorem 1.1 is a large data resulty; indeed we make no smallness assumptions 
on the constants A and TZo- 

Our theorem is connected and partially motivated by the following three earlier 
breakdown criteria results: 

1. The first is a result of M. Andersson, [And], who showed that a breakdown can 
be tied to the condition that 



Our result can be viewed as complimentary. It is clear however that the condition 
(4) is formally weaker as it refers only to the second fundamental form k and the 
lapse n which requires one degree less of differentiability Moreover a condition 
on the boundcdness of the L°° norm of R exhausts all the dynamical degrees of 
freedom of the equations. Indeed, once we know that ||R(£)||l°° is finite, one can 
find bounds for n, Vn and k on E t purely by elliptic estimates. This is certainly 
not true in our case. 

2. Our result can be also compared to the well known Beale-Kato-Majda, [BKM], 
criterion for breakdown of solutions of the incompressible Eulcr equation 

d t v + (v ■ V)u = — \7p, div v = 0, 

with smooth initial data at t — t - A routine application of the energy estimates 
shows that solution v blows up if and only if 



J t 

The Beale-Kato-Majda improves the blow up criterion by replacing it with the 
following condition on the vorticity w = curl v: 



||R(t )|| L2(Eo) <K 



(8) 



lim sup ||R(i)|| L oo = oo. 

t->tz 




(9) 




to 



(10) 



4 



SERGIU KLAINERMAN AND IGOR RODNIANSKI 



To relate Vi> and lu one observes that 

div v — 0, curl v = lu 

forms an elliptic system for v in terms of lu. Thus Vw can be expressed in terms lu 
via a singular integral operator, i.e. a zero order pseudodifferential operator: 

Vv = P°(lu). (11) 

Although P° does not define a bounded map L°° — ► i°° it can be shown that (11) 
is sufficient to reduce the breakdown condition (9) to the more satisfying one (10), 
in terms of the vorticity alone. 

Similarly, in the case of the Einstein equations energy estimates, expressed relative 
to a special system of coordinates ( such as wave coordinates) , show that breakdown 
does not occur unless 




||0g(i)|| L ~dt = oo. 



This condition however is not geometric as it depends on the choice of a full coor- 
dinate system. Observe that both the spatial derivatives of the lapse Vn and the 
components of the second fundamental form, fc^ = — in -1 dtg%j, can be viewed as 
components of dg. 

Note however that after prescribing k and Vn we are still left with many more 
degrees of freedom in determining dg. The fundamental difficulty that one needs 
to overcome is that of deriving bounds for R using only bounds for || V log n(t) ||z,oo + 
||fc(t)||i,oo and geometric informations on the initial hypersurface So- Clearly this 
cannot be done by elliptic estimates alone. Thus, as opposed to both the results of 
M. Anderson and Beale-Kato-Majda, it is far less obvious that a condition such as 
(4) can cover all dynamic degrees of freedom of the Einstein equations. Despite the 
formal similarity with the previous results mentioned above, the proof of Theorem 
1.1 requires a conceptually different treatment. 

3. Finally, the result whose proof is closest in spirit to ours and which has played 
the main motivating role in developing our approach, is the proof of global regularity 
of solutions of the Yang-Mills equations in R 3+1 by Eardley and Moncrief, see [EMI], 
[EM2] . To explain the connection of their result to ours we review below its main 
ideas. 

Recall that the curvature tensor F a p of a Yang Mills connection A a dx a , with values 
in the Lie algebra su(N) is a critical point of the Yang-Mills functional 

YM[F] = I Tr(*FAF) 

and verifies the wave equation, 

□ (A) F = F*F, (12) 
where D^ X ) denotes the covariant wave operator, 

□ (A) F = D a D a F = DF + [A, OF] + [d\, F] + [A, [A, F]], 
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□ denotes the usual D'Alembertian in R 3+1 and D a = d a + [X a , •] the gauge co- 
variant derivative. Since the Minkowski space-time metric 

m = — dt 2 + 5ijdx l dx^ 

is static (in particular n — 1 and k = 0) the energy of F associated with the 
energy-momentum tensor Q[F] Q/ g = F^F^a + *F Q A *F^a and vectorficld T = d t is 
conserved. In particular, the flux of energy T v through the null boundary N~{p) 
of the domain of dependence J~ (p) of an arbitrary point p can be bounded by the 
energy of the initial data which we denote by lo- We assume that smooth data for 
F is prescribed at t = and restrict J~(p) and Af~(p) to t > 0. We recall that 

the flux has the form, T v = ( J M - {p) Q[F](L, T) ) 1/2 with L = —d t + d r the null 
geodesic generator of JV~(p) normalized by the condition < L, T >= 1. 

The proof of the global regularity of solutions of the Yang-Mills equations is based 
on the boundedncss of the flux T p < Iq < oo. Here are a summary of the main 
steps. 

1. Rewrite (12) in the form DF = F * F - (D (A) - D)F. Using the explicit rep- 
resentation, in M 3+1 , of solutions to the inhomogeneous wave equation, we deduce, 
for all points p, with t > 0, 

F(p) = (47T)- 1 f r^FiF + F^^i) (13) 

JM- (p:5) 

- (^r 1 f r- 1 (D (A) -D)F. 

Here N~ (p, 5) represents the portion of the null cone J\f~ (p) included in the time 
slab [t(p) — 5,t(p)), with t(p) the value of the time parameter at p. Also r is the 
distance, in cuclidean sense, to the vertex p and F^ ) (p; S) represents an homoge- 
neous solution to the wave equation whose initial data at t = t{p) — S coincide with 
those of F. 

2. Ignore, for a moment, the presence of the third term on the right hand side of 
(13). Using the explicit form of the nonlinear term F*F one notices that at least 
one component of the product can be estimated by the flux T p of F through the 
null hypersurface N~ (p) . Denoting |F| = J2 a p \^ct@\, we have by a simple estimate, 

\(F(p)-F^(p;S)\ < T p ([ r- 2 ) 1/2 ||F|| L » ( ^- (p;4)) 

JM- (p;8) 

< 5 1/2 ^ 7 p ||F|| I ,oo( i7 -( Pi5 )) 

where IIFH^^-^^) denotes the sup- norm of |F| for all points in the domain 
of dependence J~{p) of p intersected with the slab t{p) — 5,t(p)]. Therefore, we 
deduce the following, 

Lemma 1.2. // 5 1 / 2 ■ T v is sufficiently small, then for any t > 

^ l|F(f " <*)lk~ + l|DF(t - S)\\ L oo (14) 
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3. Arguing recursively and using the standard local existence theorem for the 
Yang-Mills system 1 , one can find bounds for all components of the curvature tensor 2 
F(p) depending only on the fact that T p is uniformly bounded and the initial data 
data F(0) is smooth. 

4. One can show that (14) remains true even as we take into consideration the 
presence of the third term in (13). Consider for example what could be, potentially, 
the most dangerous term, 

f r _1 A-0F. 

Here we have to hope that we can integrate by parts to transfer the derivative 
from F to A. This can only be done if A • d is tangential to the light cone Af~(p). 
Miraculously, this can be achieved by taking A in the Cronstrom gauge, that is one 
assumes that the conection 1-form A satisfies 

(x - y) a \ a = 0, (15) 

where x a are the space-time coordinates of p and y a those of a point q G M~ (p) . 
With this choice, after the integration by parts, one can treat all the remaining 
terms in (D(a) — 0)F roughly in the same way as the main term F*F. To show this 
one has to observe that the value of A at any point q in the domain of dependence 
of p can be estimated by ||F||i,«>( > 7-(p,5))- This leads to the same estimate (14) as 
stated in the Lemma above. 



5. In [Kl-Ma] the global regularity result was reproved by strenghening the clas- 
sical local existence result to A G H 1 (M 3 ) and E € i 2 (R 3 ), which is at the same 
regularity level as the energy norm. That required, instead of the pointwise esti- 
mates (14), a new generation of L A type estimates, called bilinear. The premise of 
the [Kl-Ma] approach was the fact that, once we have a local existence result which 
depends only on the energy norm of the initial data, global existence can be easily 
derived by a simple continuation argument. 

6. In [Kl-Ro5] we have developed a gauge independent approach to the proof of the 
Eardlcy-Moncrief result. The approach is based on a Kirchoff-Sobolev parametrix 
for D(A)i similar to the one we use in this paper, which replaces (13) by a gauge 
invariant formula depending, implicitly 3 , only on the values of F along Af~(p). 

This paper was motivated in part by the desire to adapt the Eardley-Moncrief 
argument 4 to General Relativity. The above discussion indicates that the Eardley- 
Moncrief proof relies on two independent ingredients: conservation of energy and 
pointwise bounds on curvature, which depend only on the flux and initial data. 
Since the analogue of the Yang-Mills energy in General Relativity (the Bel- Robinson 
energy) is not conserved one can only hope to reproduce the second part of the 



In a given gauge such as the Coulomb gauge 
2 Once bounds are established for F(p) one can proceed in the same manner to derive bounds 
for derivatives of F at p. 

^Through transport equations along the null boundary of the causal past of p 
Adapting [Kl-Ma] to General Relativity is the goal of the bounded L 2 -curvature conjecture, 
see [Kl]. 
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Eardley-Moncrief argument and prove a conditional regularity result which states, 
roughly, that smooth solutions of the Einstein equations, in vaccum, remain smooth, 
and can therefore be continued, as long as an integral quantity, we call the flux of 
curvature, remains bounded. A possibility of such a result first became apparent 
to us in a discussion with V. Moncrief 5 . Such a result could also be deduced, in 
principle, from the stronger bounded L 2 -curvature conjecture, accoding to which 
the initial value problem is well posed for initial data sets with L 2 bounds on 
its curvature. In this paper we actually take a step closer to implementing the full 
analogue of the Eardley-Moncrief result. Rather then imposing a direct condition on 
the finiteness of the Bel- Robinson energy and curvature flux we formulate conditions 
(perhaps more natural albeit more restrictive) which control the extent to which 
the energy is not conserved. These conditions, which form our breakdown criterion, 
involve uniform bounds on the second fundamental form k and derivatives of the 
lapse n. 

In what follows we give a short summary of how the mains ideas in the proof of the 
Eardley-Moncrief result for Yang-Mills can be adapted to GR. 

1. A The curvature tensor R of a 3 + 1 dimensional vacum spacetime (M, g), see 
(1), verifies a wave equation of the form, 

□ g R = R*R (16) 

where D g denotes the covariant wave operator D g = D"D a . 

2. The Bel- Robinson energy- momentum tensor has the form 

Q[R]a/37<5 = RaA7^R ) 3 $ + ^a\~fn R^ §■ 

and verifies, D <5 Q Q( 3 7< 5 = 0. It can thus be used to derive energy and flux estimates 
for thee curvature tensor R. As opposed to the case of the Yang-Mills theory, 
however, in General Relativity the background metric is a dynamic variable itself 
and thus does not admit, in general, Killing fields (and in particular a time-like 
Killing field). This means that we can not associate conserved quantities to a 
divergence free Bel-Robinson tensor. It is at this point where we need crucially 
our approximate Killing condition A2. Indeed that condition suffices to derive 
bounds for both energy and flux associated to the curvature tensor R. Using the 
Bel-Robinson energy momentum tensor Q the energy associated to a slice S t is 
defined by the integral 

£{t) = f Q[R](T,T,T,T) (17) 

Js t2 

while the flux, through the null boundary M~ (p) of the domain of dependence (or 
causal past) J~ (p) of a point p, is given by the integral 

T-(p) = (j Q[R](L,T,T,T))' (18) 
JN-(v) 



5 V. Moncrief has been independently pursuing the analogy between the Einstein and Yang-Mills 
equations by developing an integral representation of the curvature tensor in General Relativity 
based on the Hadamard-Friedlander method (as in [Fried]), see [M]. 
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where L is the null geodesic generator of Af (p) normalized at the vertex p by 
<L,T>=1. 

As in the case of the Yang-Mills equations it is precisely the boundedness of the 
flux of curvature that plays a crucial role in our analysis. In General Relativity 
the flux takes on even more fundamental role as it is also needed to control the 
geometry of the very object it is defined on, i.e. the boundary of the causal past of 
p. This boundary, unlike in the case of Minkowski space, is not determined a-priori 
but depends in fact on the space-time we are trying to control. 

3. In the construction of a parametrix for (16) we cannot, in any meaningful way, 
approximate D g by the flat D'Alembertian □. To deduce a formula analogous to 
(14) one might try to proceed by the geometries optics construction of parametrices 
for D g , as developed in [Fried]. Such an approach would require additional bounds 
on the background geometry, determined by the metric g, incompatible with the 
limited assumption A2 and the implied finiteness of the curvature flux. We rely 
instead on a geometric version, which we develop in [Kl-Ro5], of the Kirchoff- 
Sobolev formula, in the spirit of that used by Sobolev in [Sob] and Y. Choquet- 
Bruhat in [Br] 6 . Applying that formula to equation (16) we obtain the following 
analogue of the formula (14): 

R(p) = - [ A-(R*R)+£+/ Err-R (19) 

where A is a 4-covariant 4-contravariant tensor defined as a solution of a transport 
equation along N~ (p, S) with appropriate (blowing-up) initial data at the vertex 
p, N~{p;5) denotes the portion of the null boundary M~{p) in the time interval 
[t(p) — S, t(p)] and the error term Err depends only on the extrinsic geometry of 
M~{p; 6). The term £ depends, in principle, only on the properties of the space-time 
in the interval [t(p) — 5,t(p) — 5/2]. 

4. As in the Yang-Mills setting the structure of the term R * R allows us to 
estimate one of the curvature terms by the flux of curvature: 

\f A-(R*R)| < r~{p) ■ \\R\\ L oo {N - M) ■ \\A\\ L2 ^- M) (20) 
JN-(p,5) 

< 5 1 ' 2 -T-{p) ■ ||R|Uoo (A f- (p , 5)) , 

provided that, 

\\M\l^-( P ;6))<S 1/2 . (21) 

Negelecting, for a moment, the third integral in (19) we can thus expect to prove a 
result analogous to that of Lemma 1.2, see proposition 5.11. 

Theorem 1.3. There exists a sufficiently small 6 > and a large constant C , 
depending only on Aq in assumptions Al and A2 as well as TZq in (8) such that 



°It is extremely important that the error term generated by our parametrix depends only on 
the geometry of the boundary of the causal past of a point. This feature is absent in all previous 
constructions. 
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for all to < t < t* , 

IIRWIU-S^C 



sup 

t-2S<t'<t-5/2 



(||R(t')||L 2 + ||DR(t)|| L2 + ||D 2 R(t)|| i2 ) 



(22) 



5. The proof of (20) depends on verifying (21). In addition, to estimate the third 
term in (19), we need to provide estimates for tangential derivatives of A and other 
geometric quantities asssociated to the null hypersurfaces Af~(p). In particular, it 
requires showing that VV~ (p) remains a smooth (not merely Lipschitz) hypersurface 
in the time slab (t(p) — S, t(p)] for some 5 > dependent only on the constants A 
and IZo. Thus to prove the desired theorem we have to show that all geometric 
quantities, arising in the parametrix construction, can be estimated only in terms 
of the flux of the curvature T~ along J\f~{p) and our main assumption Al. Yet, 
to start with, it is not even clear that we can provide a lower bound for the radius 
of injectivity of N~(p). In other words the congruence of null geodesies, initiating 
at p 1 may not be controllable 7 only in terms of the curvature flux. Typicaly, in 
fact, lower bounds for the radius of conjugacy of a null hypersurface in a Lorentzian 
manifold are only available in terms of the sup-norm of the curvature tensor R along 
the hypersurface, while the problem of short, intersecting, null geodesies appears 
not to be fully understood even in that context. The situation is similar to that 
in Riemannian geometry, exemplified by the Cheeger's theorem, where pointwise 
bounds on sectional curvature are sufficient to control the radius of conjugacy but 
to prevent the occurrence of short geodesic loops one needs to assume in addition 
an upper bound on the diameter and a lower bound on the volume of the manifold. 

In a sequence of papers, [Kl-Rol]-[Kl-Ro3], see also [Wang] 8 we have proved lower 
bounds on the geodesic radius of conjugacy of null hypersurfaces. The methods 
developed in those papers can be adapted to also prove lower bounds on the radius 
of conjugacy with respect to the time parameter 9 t. It may however be possible 
that the radius of conjugacy of the null congruence is bounded from below and yet 
there are past null geodesies form a point p intersecting again at points arbitrarily 
close with respect to the time parameter t) , to p. In [Kl-Ro4] we have shown that 
this cannot happen in a space -time verifying our conditions Al and A2. Thus 
the combined results of [Kl-Rol]-[Kl-Ro4] allow us to derive a lower bound on the 
radius of injectivity of J\f~{p) depending only on A . 

6. As in the case of Yang-Mills equations one can use the result of Lemma (1.3), 
together with the classical local existence result for the Einstein equations, such as 
that in [C-K], to show that solutions can be extended as long the bounds on ( T '7r 



Different null geodesies of the congruence may intersect, or the congruence itself may have 
conjugate points, arbitrarily close to p. 

8 In [Kl-Rol]-[Kl-Ro3] we have considered the case of the congruence of outgoing future null 
geodesies initiating on a 2-surfacc So embedded in a space-like hypersurface So- The extension 
of our results to null cones from a point forms the subject of Qian Wang's Princeton 2006 PhD 
thesis, see [Wang]. 

9 The results in [Kl-Rol]-[Kl-Ro3] and [Wang] were proved with respect to the geodesic foliation. 
In this paper, as well as in [Kl-Ro4], we rely on an extension of these results to the foliation on 
Af~(p) induced by the space-like foliation St. 



hold true. 
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Finally we would like to point out possible refinements of our main theorem 1.1. 
We expect that one should be able to replace the pointwise condition A2 with the 
integral condition, 

f ||7r(t)||ioodt < oo (23) 

J to 

Moreover it may be possible to improve the result even further by eliminating the 
term Vlogn in (4) or (23) and requiring instead only a pointwise bound on n. 



2. Constant Mean Curvature foliations 

As described in the introduction (A4*, g) denotes a Lorentzian manifold of the form 
M.* = IxE, where E is a three dimensional, compact, connected, orientable smooth 
manifold foliated by a CMC foliation £ t with lapse n and second fundamental k, 

n= ( - g(T>t,T>t)y 1/2 , k(X,Y) = g(D x T,Y) 

where T denotes the future unit normal to S t . The time interval I = [to, t*), where 
t = —1 and t* < 0. 

We decompose a space-time vcctorfield X relative to the unit timelike T, 

X = X a T + X, <T,A>=0, (24) 

We define the positive definite Ricmannian metric, 

h(X,Y)=X°-Y°+g(X,Y). (25) 

where g denotes the metric induced on X t . We can also write (25) in the form, 

h a ^g aP + 2T a T l3 . (26) 

Given a space-time tensor U we denote by \U\ its norm with respect to the metric 
h. More precisely, if U is a to— covariant tensor, 

\U\ 2 = h iljl ■ ■ ■ h imim U il ... im U h ...j m (27) 

The following bound follows immediately from our main assumption (5), 

|DT| <A (28) 

Since D 7 /i a/3 = 2(D 7 T a T /3 + T a D 7 T /3 ) we have \Dh\ < 4|DT|. Therefore, 

\Dh\ < A (29) 

Also, since the components of the deformation tensor ir = ^tt = CtS are given 

by 

7T o = 0, 7t 0j = n _1 V 4 n, ir i3 ■— n~ x d t gi 2 ■ — -2k i3 ■, 

we have, 

l (T M<A (30) 

Given two tensors U, V we shall denote by U ■ V any tensor which is obtained from 
the tensor product of U and V by taking contractions with respect to the space-time 
metric g. Clearly, 

\u ■ v\ < \u\ ■ m 
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For any coordinate chart O, with coordinates x = (x , x 2 , x 3 ), we denote by (x = 
t, x 1 , x 2 , x 3 ) the transported coordinates on/xO obtained by following the integral 
curves of T. In these coordinates the metric g takes the form 

g = -n 2 dt 2 + g lj dx t dx J , (31) 

Relative to these coordinates t, x we have the equations, 

dt9i = -2nk tj (32) 
d t kij = —ViVjn + n(Rij+tYgkkij—2k ia k a j) (33) 

with Rij the Ricci curvature of of the induced metric g on E t . We also have the 
constraint equations, 

R - \k\ 2 + (trk) 2 = (34) 

V J '% = Vitrfc. (35) 

In view of the constant mean curvature condition on the foliation E t we can always 
reparametrize t so that 

tr g fc = t. (36) 

As mentioned in the introduction we can assume that the initial hypersurface E 
corresponds to the value t = t — — 1. In view of (33), (34) and (36) we deduce the 
lapse equation, 

An=\k\ 2 n-l (37) 

At a point p of minimum for n we must have \k\ 2 n — 1 > 0. Therefore, at p 
n > \k{p)\~ 2 . On the other hand, since \k\ 2 = |fc| 2 + i(trfc) 2 , at a point of maximum 
we have, \k\ 2 n + |(trfc) 2 n — 1 < 0. Therefore, 

< n < i. (38) 



\T, t \ ~ — ( \J det gdx = — f nty/det gdx. 



ll*(*)llicc - ~t 2 - 

Observe also that, since d t log(det g) = — 2ntrfc = — 2nt, 
d_ 
dt 

where |E t | denotes the volume of the compact manifold E t . Thus, 

o < ^li:*! < ai^ntr 1 

As a consequence of (36) the ratio of the volumes of |S t | and |S | can be estimated 

by, 

i < M < l*oi 3 



Sol - |t| 3 

Therefore, since to = — 1, we have proved, 

Proposition 2.1. For all — 1 = to < t < t* < we have the bounds, 

mm:- 71 -"* (39) 
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Moreover, i/|£t| denotes the volume of St and £o = £t , 

|So| < |S t | < i^lEol < \t\~ 3 (40) 

2.2. Coordinate estimates. We recall the following lemma, see lemma 2.2 in 
[Kl-Ro4]. 

Lemma 2.3. I/E * s compact and verifies Al of the introduction, there must exist 
a number p a > such that every point y E E admits a neighborhood B, included 
in a neighborhood chart U, such that B is precisely the Euclidean ball B = Bpf (y) 
relative to the local coordinates in U . 



Proof : See the proof of Lemma 2.2. in [Kl-Ro4]. ■ 
Next we recall the result of proposition 4.1. in [Kl-Ro4]. 

Proposition 2.4. If assumptions Al and A2 are verified, then there exists a large 
constant C = C(A ) such that, 

C-'ltf <g i3 (t,x)Ce <C\e, Vxe(7 (41) 

Proof : For convenience we reproduce the proof given in [Kl-Ro4]. We fix a 
coordinate chart U and consider the transported coordinates t, x 1 , x 2 , x 3 on I x U. 
Thus d t gij — —\nkij. Let X — X be a time-independent vector on M tangent to 
S t . Then, 

d t g{X,X) = - l -nk{X,X). 

Clearly, 

\nk{X,X)\<\nk\ g \X\l<\\nk{t)\\ L ~\X\l 

with \k\ 2 g = g ac g bd k ab k cd and \X\ 2 g = .Y'.Y % = g(X, X). Therefore, since d t \X\ 2 g = 
dtg(X,X), 

- l -\\nk{t)\\ L ~\X\ 2 g < d t \X\ 2 g < \\\nk{t)\\ L ~\X\l 

Thus, 

\X\ ga e-^ W nk ^™ dT < \ X \ 2 gt < \X\ go e f *o ^)\\L~dT 
from which (41) immcdiatley follows. ■ 



2.5. Sobolev inequalities. The properties of local transported coordinates estab- 
lished in the previous section can be used to prove the following Sobolev inequality 
for scalar functions. 

Proposition 2.6. Assume assumptions Al and A2 verified. There exists a con- 
stant C depending only on Aq such for every smooth scalar function on E t , to < 
t <t* such that 

ll.fllL3( St) <C(||V/|| L1(St) + ||/|U 1(Et) ) (42) 
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Proof : By a partition of unity we may assume that / has compact support in a 
local chart V = St n (/ x U) of transported coordinates t, x = (a; 1 , a; 2 , a; 3 ). Then, 
writing 

/X 1 rx 2 px £ 

dif(y,x 2 1 x 3 )dy = / d 2 f(x 1 1 y,x 3 )dy = / d 1 f{x 1 ,x 2 1 y)dy, 
-OO J — OO J — OO 

\f(x)\ 3/2 < ( [ X \&if{y,^,^)\dy f \d 2 f{x\y,x 3 )\dy f Id^ix^x^y^dy) 

\J—oo J —oo J —GO / 

Thus, by Holder, 

/ \f(x)\ 3/2 dx<([ \Vf(x)\dxf/ 2 
Jv Jv 

Therefore, since in view of (41) we have C" 3 / 2 < J\g\ < C 3 / 2 , 



1/2 



(/ \f{x)\ 3/2 V\g\dx)* < / |v/(x)|/M^- 

J v Jv 
which proves (42) as desired. Inequality (42) is an immediate consequence of (42). 



Corollary 2.7. For any smooth tensorfield F on St and any 2 < p < 6, 

||F||L- (Et ) < C(\\WF\\% 2 { -f\\F\\% p ( -f + \\F\\ L2{St) ) (43) 

Proof : We have, 

\\F\\l P J 3 = |||F| 2 f/ 3 || L 3 /2 <C(||V|F| 2 f/ 3 || il + |||F| 2 P/ 3 || L1 ) 

< ciwwFW^ + WFW^-wiFi 2 ^ 3 - 1 ^ 

< C(\\VF\\ L2 + \\F\\ L2 ).(\\F\\ Li£T e)^ 
In the particular case when p = 6 we derive, 

ll^llla<C(||VF|U 2 + ||^|U 2 ).||^||l 6 

Therefore, 

\\F\\ L e<C(\\VF\\ L . + \\F\\ L .) 

Similarly, for p = 3, 

||F|| 2 3 <C7(||V^|U 2 + ||F|| L2 ).||F|| i2 

and thus, 



||F|| L3 <C(||VF|| i2 + ||F|| i2 ) 1/2 .||F||^ 2 
The general case follows by interpolation. 

Here is another useful simple calculus inequality which we will make use of. 
Lemma 2.8. Let F be a tensorfield on a compact Riemanian manifold. Then, 

||VF|U, < ||V 2 F||^ 2 ||F||^ 
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Proof : After an integration by parts and Holder, 

|jVF| 4 <||V 2 f|| L2 ||V^||i 4 ||F|| L » 

Hence, 

ltVF||| 4 < \\V 2 F\\ L 2\\F\\ L oo 



Remark. We cannot use transported coordinates to derive a Sobolev inequality 
of the form, 

ll/llL»<l!v 2 /|| i2 + ||/|| L2 

even in the case of a scalar function /. Indeed, the standard Sobolev inequality in 
a coordinate chart U provides, 

3 

ll/IU-w < E II WIU'oo + Wfhnu) 

ij = l 

On the other hand VjVj/ = didjf — T^dif and therefore we cannot derive the 
desired estimate without a bound for the L 3 norm of T. Unfortunately, the only way 
to estimate T is by differentiating the equation d t g = —Ink from which we could 
only bound its L 2 norm. To get around this difficulty we need a better system of 
coordinates. In [Kl-Ro4] we have a proved a slightly more general version of the 
following: 

Theorem 2.9. Assume that A4* is globally hyperbolic and verifies the assumptions 
Al and A2 as well as (8). Then, for any e > 0, there exists ro > 0, depending 
only on e, Ao, IZo, i*, such that on any geodesic ball B r C T,t, r < ro, centered at 
a point p t £ St, there exist local coordinates relative to which the metric g t verify 
conditions 

(1 + €)-% < 9ij < (l + e)«% (44) 
r [ \d 2 g tl \ 2 dv g < e. (45) 

JB r (p) 

As a corollary we derive the following version of the Sobolev inequlity 
Corollary 2.10. Given a smaooth scalar function f on S t we have, 

ll/IU~ (Et ) < C\\V 2 f || i2(Et) + ||/|U 2(St) ) (46) 
with C > a universal constant, i.e. depending only on the fundamental constants 

3. Basic Curvature energy estimates 

3.1. General procedure. We recall the general procedure to derive energy esti- 
mates for R, see section 7.1 in [C-K]. First let W denote a Weyl field, i.e a four 



CURVATURE 



15 



covariant tensor traceless tensor W a /3js verifying all the algebraic symetries of the 
curvature tensor R. Let, 

Q[W] a07 s = W aX ^w£f + *W aXjfl *Wj£f (47) 

Given a vectorfield X we denote P a = Q[W] a /3 1 sX ,3 X 1 X s . By a straightforward 
calculation, 

T> a P a = T)WQ[W}X f) X''X s + ^Q Q/ 3 7 «5 {x) Tr af> X~<X 5 

where ( x ^n is the deformation tensor of X Therefore, integrating on the slab 
UfG[t ,t]S t / we derive the following. 

Proposition 3.2. Let Q = Q[W] be the Bel-Robinson tensor of a Weyl field W. 
Then, 

j Q(X,X,X,T) = I Q(X,X,X,T)+ f f DivQ(X, X, X)ndv g 

+ If I <W {X) 7T afj X^X s ndv g (48) 
with dv g denoting the volume element on S t . 

The following proposition is an immediate consequence, see also section 5 in [Kl-Ro4] 
for a proof. One simply needs to apply the proposition above for X = T together 
with the positivity of Q(T, T, T, T) and the uniform bounds for ( T '7r and n. 

Proposition 3.3. Under assumption A2 There exists a constant C — C(Ao,t*) 
such that, for any to < t < t* < 0, 

\\K(t)\\ L 2 <CUo. (49) 

where IZo is the constant defined by (8). 

Definition 3.4. In what follows we extend the usual notation A < B to include 
inequalities A < cB where c = c(t*, Ao, TZo) is a constant which depends on our 
fundamental constants t*, Ao and TZq. 

In particular, in view of proposition 3.3 we can write 

||R(t)|| ia <fto<l. 

3.5. Wave equation for the curvature tensor. Recall the Bianchi identitities, 

D [<T R a/3]75 = (50) 

or, equivalently since R Q) 3 = 0, 

D 5 R Q/3 ^ = (51) 
Differentiating (50) once more an taking the trace we derive, 
□ R a ^ 7 5 + D <J D ct R^ cr7 5 + D cr D^R crct7 ,5 = 
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Now, in view of (51), commuting covariant derivatives, 
Hence, 

D cr D ct R ( 3 (J7( 5 + D cr D ( 3R (JQ7 5 = R^ CT7 5 (R^ CT Q — R a MCT g) 

— R acr7M (R (5 At (3 + R^ ^) 

— Racr^(5(R ( 3 ^ 7 + R 7 ^ £j) 

Thus introducing the notation, 

(R*R) Q( 3 7 ,5 = — R M<T7 5 (R^^ — Rj/"^) + R Q (t 7M (R 5 M< ^ + R/j^) (52) 

+ Racr^i (R-^^ + R- 7 M °a) 

we derive, 

□ R = R* R (53) 

Clearly W = R * R is a Weyl field, i.e. it satisfies all the algebraic symmetries of 
the curvature tensor plus the traceless condition Wj 1 p = 0. 

3.6. Energy estimates for higher derivatives. To estimate the first derivatives 
of R we shall use the covariant wave equation (53). Recall the positive definite 
space-time metric h defined by (25). Given a tensor-field U ai ... am we write, for 
simplicity, 

h IJ U T Uj = h a ^ ...h a ^U ai ... am U (il ... f)m 

Ul = U ai ... am , Uj = Up 1 ... 0m , h" = h a ^...h a ^ 

Consider the energy-momentum type tensor Q^g associated with the covariant 
wave operator □ acting on tensors, 

Q^[UU: = h IJ T> a U I T>( 3 Uj-^g af3 h IJ gL""Dt l U I T> v Uj (54) 

We have, 

D*QM [U] aP = h'^Uj^-DpUj + h IJ D^T> a UjlDpUj - g^/i 7 V^D^D.t/, 

= h IJ B a U I (DU J ) + h IJ (D^ a U! - B^pU^Uj 

+ D^h IJ {T> a UiT> Uj - ig^g^D^C/.D^,;) 
Consequently, in view of (29), 

|DQW[C/]| < \DU\\UU\ + \R\\U\\BU\ + A \BU\ 2 
Therefore since, 

D^(Q (ffi) [C/]^r) =^T a QW[[/]^ + TWQW[[f]^ (55) 

we derive, 

|D^Q^[[/] Q/3 T Q )| < |D(7||n(7| + |R||C/||DC/| + Ao|DC/| 2 (56) 
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On the other hand, 

k 

Q W P](T,T) = ^"(Dotf/Dotfj + ^DjtfjD^j) 

1=1 

= \\vu\ 2 

Integrating (55) we derive, 

/ \DU\ 2 = 2 [ QM[[/](T,T) 
Js t Js t 

I \DU\ 2 +ff \D' 3 {q} w \U] aP T a )\ 
f f (|DC/||DC/| + |R||;7||DC/| + Ao|D[/| 2 ). (57) 
Applying this to U — R and using the equation (53) we obtain, 

||DR(t)||£ a < ||DR(io)||| 2 +A f \\T>R(t')\\l 2 dt> 

J to 

+ ADR(tOIM|R(OIM|R(*')IU^' 

Jt„ 

< ||DR(io)||£ 2 +A /"||DR(i')||| 2 d*' 

Jto 

+ TZo [ ||DR(t , )|| i a||R(t')IU~dt' 



< 



< 



Jt 

Therefore in order to get an a-priori estimate for j|DR(i)|| L 2 it suffices to prove an 
estimate for the L°° norm of R. More precisely, 

Proposition 3.7. Assume that the assumptions Al, A2 hold true. Then the 
following derivative curvature estimates hold true for all to < t < t* , 

||DR(i)||| 2 < C(||DR(t )||| a + /* ||R(i')||!oodO (58) 

Jto 

with C a constant depending only on A , TZo an d t*. 

To estimate the second derivatives of R we apply (57) to the tensor U = DR. 
Thus, 

f |D 2 R| 2 < f j (|D 2 R| |D(DR)| + |R||DR||D 2 R| + A |D 2 R| 2 ) . 
< [ f (|R||DR||D 2 R| + A |D 2 R| 2 ) 

Jto JS t , 
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Hence, 



|D 2 R(t)|| 2 < 



D 2 R(i )|| 2 2 +A / ||D 2 R(t')||!^' 

J to 

+ f 1 1 D 2 R(t') || i2 ||DR(t') || i2 1| R(t') lk~ rf*' 

J to 

< ||D 2 R(t )||| 2 +A /*||D 2 R(t' 

Jtn 



'^" 2 dt' 



+ f ||D 2 R(t')||L 2 ||DR(t')llLH|R(^)IU~*'- 

J to 

We therefore, deduce the following, 

Proposition 3.8. Under the same assumptions as in proposition 3.8 we have, 

||D 2 R(t)|| 2 2 < C(||D 2 R(t )|| 2 2 + /"* ||r>R.(f )lli a ||R.(* / )lli«dt / ) (59) 

J to 

with C a constant depending only on A , 7£o an d t*- 

4. Past null boundaries 

The goal of this section is to review the main result of [Kl-Ro4], concerning the 
null boundaries of past domains of dependence, and show how they apply to our 
situation. Starting with any point p in a subset M* — U te [ t0)t „)S t of M, we 
denote by J~{p) — 3~ (j>; M.*) the causal past of p, relative to A4*, by I~(p) its 
interior and by Af~ (p) its null boundary. In general Af~ (p) is an achronal, Lipschitz 
hypersurface, ruled by the set of past null geodesies from p. We parametrize these 
geodesies with respect to the future, unit, time-like vector T p . Then, for every 
direction w 6 § 2 , with S 2 denoting the standard sphere in R 3 , consider the null 
vector £^ in T p M, 

g(4,,T p ) = l, (60) 

and associate to it the past null geodesic 7^(s) with initial data 7^(0) = p and 
7w(0) = tv We further define a null vectorfield L on J\f~ (p) according to 

H"tu,(s)) = 7^(s). 

L may only be smooth almost everywhere on 7V~ (p) and can be multi- valued on 
a set of exceptional points. We can choose the parameter s in such a way so that 
L = 7u>(s) is geodesic and L(s) = 1. 

For a sufficiently small S > the exponential map Q = Q ~ defined by, 

(s,w) ->7o,(s) (61) 

is a diffcomorphism from (0,8) x § 2 to its image in J\f~(p). Moreover for each 
io G § 2 either 7 w (s) can be continued for all positive values of s 10 or there exists a 
value s»(u>) beyond which the points 7a, (s) are no longer on the boundary Af~(p) 
of J~{p) but rather in its interior, see [HE]. We call such points terminal points 



10 for which J w (s) stays in M„ 
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of 7V _ (p). We say that a terminal point q — 7 w (s*) is a conjugate terminal point 
if the map Q is singular at (s*,u>). A terminal point q = 7 w (s„) is said to be a cut 
locus terminal point if the map Q = Q~ is nonsingular at (s*,w) and there exists 
another null geodesic from p, passing through q. 

Thus Af~ (p) is a smooth manifold at all points except the vertex p and the terminal 
points of its past null geodesic generators. We denote by T~(p) the set of all 
terminal points and by Af~ (p) = Af~ (p) \ T~ (p) the smooth portion of Af~ (p) . 
The set £/ _1 (T~(p)) has measure zero relative to the standard measure dsda§2 
of the cone [0, oo) x § 2 . We will denote by cL4jvw p ) the corresponding measure 
on J\f~(p). Observe that the definition is not intrinsic, it depends in fact on the 
normalization condition (60). 

Definition 4.1. Given p G Al» we define i*(p) to be the supremum over all 
the values s > for which the exponential map Gp '■ (s,w) — ► 7 w (s) is a global 
diffeomorphism. We shall refer to i~(p) as the past null radius of injectivity at p 
relative to the geodesic foliation defined by (60). 

We also define i~(p,t) (the null radius of injectivity relative to the t-foliation) to 
be the supremum over all the values t(p) — t, t < t(p), for which the exponential 
map g = Q~ t , 

(t,w) -7u,(*) =7« (*(*)) (62) 
is a global diffeomorphism. 

Definition 4.2. We define d~(p,t) to be the distance, measured with respect to 
the time parameter t, from p to the past boundary of X» C M. 

The following theorem is an immediate consequence of the Main Theorem II proved 
in [Kl-Ro4]. 

Theorem 4.3. Assume that A4* is globally hyperbolic and verifies the assumptions 
Al and A2 as well as (8). There exists a positive number > 0, depending only 
on A 0; TZo, and t* < 0, such that, for all p e M*, 

i~(p,t) > min(«», cT(p,t)) (63) 

Proof : According to the Main theorem II of [Kl-Ro4] and the remark following it, 
depends only on our main constants, A , TZo and a constant which provides 
uniform bounds for the lapse n, 

A _1 < n < N Q . 

The finiteness of N , N^ 1 follows from (38) and the assumption i» < 0. ■ 

Once we have a lower bound for i*(p,t) it is straightforward to also get a lower 
bound for the radius of injectivity i*{p) with respect to the geodesic foliation. 
Indeed all we need is to show that s does not vary much (along A/"~ (p)) as a function 
of t in a time interval of size 1. This follows immediately from the following. 
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Lemma 4.4. There exists a constant c > 0, depending only on Ao, such that, 

c > < |f | < , (64) 

Proof : Wc introduce the null lapse, 

tp-^g^L) (65) 
Observe that ip > with ^s(p) = 1. Moreover 

I = -("«""' < 66 > 

with n the lapse function of the t foliation. On the other hand, we have 

L = -cp-^T + N) 
with N of length 1 perpendicular to T. Now, 

= A g(T , i)=g(DiT , i) = _i miril 



Therefore, 



from which, 



^- 2 ( (T) ^TAT + i (T) ™)). 



I^I<,- 2 A 



\<p(s) - 1| < A oS . (67) 

Thus, for an interval in s of size 1 we deduce that 2 _1 < 0(s) < 2 and therefore, in 
view of the uniform bound for n of proposition 2.1, we infer that there must exist 
a constant c > 0, depending only on A , such that (64) holds. 



4.5. Geometry of smooth null cones. In this subsection we provide additional 
geometric informations for the null boundaries N~ (p, 8) with 5 < i~(p,t) with 
K{Pi a lower bound for past null injectivity radius with respect to the t-foliation. 
Here Af~(p, 5) denotes the portion of N~{p) for t between t(p) and t(p) — S. 

Let S t denote the 2 dimensional space-like surfaces of intersection between S t and 
N~{p). At any point of JV~(p, S) \ {p} we can define a conjugate null vector L 
with g(L,L) = —2 and such that L is orthogonal to the leafs St- In addition we 
can choose {e a ) a =i,2 tangent S t such that together with L and L we obtain a null 
frame, 



g(L,L) = -2, g(L,L) = g(L,L) = 0, 
g(i, e a ) = g(L, e a ) = 0, g(e a , e b ) = S ab . 



(68) 
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We denote by 7 the restriction of g to St i.e. j(X,Y) = g(X,Y) Endowed with 
this metric St is a 2 dimensional compact riemannian manifold. We denote by y 
the restriction of D to S t , Clearly, for all X, Y e T(S t ), 

f x Y = T> X Y + ±g(B x Y,L)L + ^g(B x Y, L)L 

We recall, see [Kl-Rol] the definitions of the following basic geometric quantities: 

Definition 4.6. The null second fundamental forms x, X, torsion £ and the Ricci 
coefficient r\ of the foliation S t are defined as follows: 

Xab = g(D a L , e b ), x ab = s(D a L , e b ), (69) 
(a = ls(D a L , L), % = ig(e a , D^L). (70) 

In addition we define tr* = 1 ab Xab, Xab = Xab - \^Xlab and 

w = ~g{D L L, L), n = L(trx) + ^tr^trx + 2wtrx 

We note that 

X ab = -f 2 Xab + 2fk ab , 
V a = -(a - (ntp)- 1 e a (rup) , 
uo = tpn~ 1 N(n). 

Our conventions imply that 

f x Y = T> x Y- l -x{X,Y)L- l - X {X,Y)L (71) 

We extend the definition of y to any covariant S— tangent tensor 7r by the usual 
formula, 

f x Tt{Yi,...,Y k ) = X*{Y x ,...,Y k )-*{y x Y u ...,Y k )-...-*(Yi,...,y x Y k ) 

with X, Yi, . . . Y k 5- tangent. Given an S'-tangent vector-field X we define f ' L X to 
be the projection to S t of T>lX, 

y L x = n L x + ± g (n L x,L)L 

We extend the definition to any covariant 5-tangent tensor it by 
f L n(Y u ...,Y k ) = Ltt(Y u ... ,Y k )-n(y L Y 1 ,...,Y k )-...-ir(Y 1 ,...,y L Y k ) 
with Yi, . . . Y k S'-tangent. 

Given an S-tangent tensor 7r wewrite y-K — (yn, y L n) and 

IfH 2 = \y L n\ 2 + \yn\ 2 . 
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4.7. Tangential covariant derivatives of space-time tensors. In this section 
we make sense of covariant derivatives of space-time tensors, not necessarily S- 
tangent along a fixed surface S — S t C N~ (p, 5) . 

We start by denning a covariant derivative for space-time vector defined on S. 
Thus we view A as a section of a vector bundle T*M over S. We interprete the 
covarant derivative ^7A of A along S as a 1-form on S with values in T*M. Thus, 
for every vectorficld X e TS and any vectorficld Z in TM, 

fA(X; Z) = f x A(Z) = X(A(Z)) - A(D X Z) = B X A(Z) 

We also write, 

(f x A) M = X a T> a A^ VX e TS. 
We define )P 2 A, the second covariant derivatives of A along S, by the formula, 
f 2 A(X, Y; Z) = (f x f A)(Y; Z) = X(fA(Y; Zj) - fA(f x Y; Z) - fA(Y; D X Z) 
or, for simplicity, 

f 2 A^(X,Y) = (f Y (f x A)), (f fYX A), 

These definitions can be easily extended to higher covariant derivatives along S and 
to higher order tensors A. 

Given aAan S'-tangent 1-form on M with values in TM we define 

f h A(X-Y) = L(A(X;Y))-A(f L X;Y)-A(X;D L Y), \/X G TS, Y e TM. 

This defintion extends naturally to higher order tensors A. Note that for a scalar 
function A on M we have 

f L A = T> L A 

4.8. Commutation formula. In what follows we will need the following commu- 
tation lemma, see [Kl-Rol]. 

Lemma 4.9. Let A^ be a function on M with values in TM verifying the equation 

f L A = F (72) 
for some TM valued function F. Then, 

yj^cA) + Xa b f b A^ = f a F^ + (Ca + vJF^ + R„\ a A x . (73) 

4.10. Curvature flux. 

Definition 4.11. The curvature flux along 11 Af~{p), is defined as follows. 

r rt(p) r 

T( P )= Q[R](T,T,T,£) = / ndt Q[R](T, T, T, L)dA t 

JM-{p) Jto J S t 

with dA t the area element of St- We also let 

T(p,5)= Q[R](T,T,T,L)= / ndt Q[R](T, T, T, L)dA t 

JM-{p,5) Jt(p)-5 J S t 

1:L Givcn a scalar function / on J\f~(p) we denote its integral on M~{p) to be, Sj^-i p ^ f = 
ItT ndt 1st f dA * = f*- (p) / dA M- (P) ■ ■ 
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to be the curvature flux along Af (p, 5) for 6 < i f (p, t). 

The following is an immediate consequence of the energy estimates of section 3.1, 
see also [Kl-Ro4]. 

Proposition 4.12. Under assumptions A2 as well as (8) the flux of curvature 
Af~{p) (denoted T{p)), can be bounded by a uniform constant independent of p. 
More precisely, for all p with t < t(p) < t* < 

Hp) < c(u,a )k . 

where C is the constant of proposition 3.3 

Proof See section 5 in [Kl-Ro4]. ■ 

We can allso introduce the reduced flux, or geodesic curvature flux, 

K(p) = { / H 2 + |/?| 2 + H 2 + H 2 + |^| 2 ) 1/2 (74) 

Jt J s t 

as well as 

K(p,5) = ( / H 2 + |/?| 2 + |p| 2 + |a| 2 + |^| 2 ) 1/2 (75) 

Jt(p)-6 J S t 

where a, f3, p, a, [3 are the null components of the Riemann curvature tensor relative 
to the St foliation: 

Otab = R-LaLb , Pa = -jR-aLLL , P = ^R-LLLL 

°~ = -^*R-LLLL, §_ a = -^RaLLL, O ab = KiaLb (76) 

Proposition 4.13. Under the same assumptions as in proposition 4-12 we have, 
with a constant C depending only on t*, Aq, 

K(p) < C(t»,A )K a (77) 



Proof We can express L,L in the form, 

L = (£ _1 (T + N), L = <p(T-N) 
where N is unit normal of S t on £ t and if the null lapse defined by 65. Also, 

T = ipL + ip- x L (78) 

Therefore, 

Q(T, T, T, L) = ^Q(L,L,L,L) + 3ip 2 Q(L,L,L,L) 
+ 3ipQ(L,L,L,L) + Q(L,L,L,L) 
= ^H 2 +3^ 2 |/?| 2 + 3^ 2 + M 2 ) + |£| 2 

and the result follows from the bound (67) for <p. ■ 



24 



SERGIU KLAINERMAN AND IGOR RODNIANSKI 



We can also get additional estimates for the flux associated to the first derivatives 
of the curvature tensor. To see that we go back to the derivation of theorem 3.7. 
We now integrate (55) in J~{p) and derive, 

/ Q W [R](T,L)< f Q(™>[R](T,L)+ / |D"(Q(™)[R] Q ^T Q )|. 

Similarly, 

/ Q W [R](T,L)< f QW[R](T,L)+/ |D^(Q^[R] q/3 T q )|. 

Here J~{p) is the causal past of p and J~{p, 8 the portion of J~{p) to the future 
of £i(p)_<5. Now, 

Q W [R](T,L) = QWIRK^ + ^'I,!) 

and, 

QW[R](L,L) = |^R| 2 
Q(-)[R](L,i) = |^R| 2 
We introduce the flux quantities, 

T^(p) = [ (\fK\ 2 + \f L R\ 2 ), 

JN-(p) 

TU(p,6) = f (|yR| 2 + |y L R| 2 ) (79) 
We can therefore reformulate proposition 3.7 as follows, 

Theorem 4.14. Assume that Al, A2 hold true. Then, for any < S < i„, with 
i* > defined by theorem 4-3, 

||DR(t)|| L2 +sup.F< 1 >(p,J) < cf||DR(t-J)|| L2 + ( /* HR^Olkocrft') 1 
pes t V Jt-s 

(80) 

w^/j C a constant depending only on A , TZq, t*. 



4.15. Estimates for the Ricci coefficients. In this section we state without 
proof a proposition concerning the regularity properties of the Ricci coefficients 
t r X: Xj C an d V as well as mass aspect function /i associated to the St foliation. 
A similar result was proved for the corresponding quantities associated with the 
geodesic foliation in [Kl-Rol], see also [Wang]. The methods used for the geodesic 
foliations can be easily adapted to prove the result below. 

Proposition 4.16. For any t e (t(p)—6,t(p)) withS < i~(p,t) the Ricci coefficients 
tTX.> X> Cj V an d M satisfy the following estimates. 

sup|*rx--£r| + || sup (t(p)-t)\ftrx\\\ L 2<C, (81) 
sup/ (|x| 2 + |C| 2 + H 2 )M dt<C, U\\l* { m- M )<C 
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with a constant C depending only on Ao,£* and curvature flux 1Z(p,6)). Here 
the point Af~(p,S) are parametrized by the coordinates (t,u>) with w 6 S 2 . The 
volume forms dA^j--t p ^\ on M^{p 1 8) and dSt on St are respectively equivalent to 
the expressions (t(p) —t) 2 dt da§2 and (t(p) —t) 2 da§2 withda^ denoting the standard 
volume form on S 2 . Notation L 2 above refers to the L 2 norm with respect to the 
measure da§2. Finally, the quantities (t(p) — t) and s(t) are equivalent. 

5. KlRCHOFF-SOBOLEV PARAMETRIX 

Earlier in this paper, see Propositions 3.7, 3.8, we were able to derive L 2 estimates 
for derivatives of the curvature tensor which depend on the additional assumption 
on the boundcdness of the L°° norm of the curvature tensor. To estimate the 
latter we rely on a special version of the Kirchoff-Sobolev parametrix introduced in 
[Kl-Ro5]. 

5.1. Optical function. To make sense of our Kirchoff-Sobolev formula we need 
to define an optical function 12 u, in a neighborhood of Af~ (p, 8), < S < i*(p, t), 
such that it vanishes identically on Af~(p, 6). Here p is an arbitrary point of A4* = 
U te [j 0;t)i )£ t . We recall that we have assumed that M* is globally hyperbolic with 
Cauchy hypersurface S to . We define u uniquely relative to the time-like vector T p 
as follows: 

Let e > a small number and T e : (1 — e, 1 + e) — > M* denote the timelike geodesic 
from p such that r e (l) = p and 1^(1) = T p . From every point q of L e let Af~(q) 
be the boundary of the past set of q and JV~(q, S) defined as before. 

We now define u to be the function, constant on each JV~(q,5), such that for 
q = T(t), 

u \^-( q ) =*-!■ 

This defines a smooth function u which vanishes on N~ (p, 5) and verifies the eikonal 
equation, 

g a0 d a ud u = 0. (83) 

Observe that the vectorficld L — g a Pdpud a is null, geodesic and verifies the nor- 
malization condition, 

g(L,T p ) =T P («) = 1. 
Thus L is the same as the vectorficld L defined earlier in section 4. 

5.2. Main representation formula. We shall next state a result which was 
proved in [Kl-Ro5], concernig tensorial wave equations of the form D^f = F, with 

a fc-covariant tensorfield. Let p G M.* and S < i~(p,t). Let A be a tensor-field 
of the same order verifying, 

D L A + ^Atr X = 0, sA(p) = J on Af~ (p, S) (84) 



'i.e. a function which verifies (83) below. 
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where Jo is a fixed k-tcnsor at p, | Jo (p) | < 1- 

Theorem 5.3. Let p G A4* and S < i~(p,t). Let be a k covariant tensor-field 
vanishing identically for t < t(p) — S. Then, given A a solution to the transport 
equations (84) we have 13 , 

^(p)-3 = - f (AD*- ±A-R(-, +1 [ 

JM-(p,5) v 1 II JAT-( P ,6) 

+ I (M + Ca^A)-* (85) 



We apply the theorem to the tensor-field $ = /(i)R where < / < 1 is a smooth 
function supported in the interval [t(p),t(p) — 6} and identically equal to f in the 
interval [t(p),t(p) — <5/2]. Since R verifies (53) we have, 

□ (/R) = /R ★ R + (D/)R + 2D a /D Q R 

with R * R defined in (52). In view of the theorem above we have the formula, 

47rR(p).J = I(p) + J(p)+K(p) + L(p)+£ (86) 

J(p) = f A-/(R*R) 

JM-{p,S) 

1 



J(p) = -- / A-R(.,.,L,L)./R 

1 JM-{p,S) 

k{p) = [ (M + CafA) ■ /h 

L(P) = \ t MA-/R 

)M-{p,5) 



2 



£ = f (□/(A.R) + 2D Q /(A.D Q R)) (87) 

JM-(p,S) 



5.4. Estimates for /(p). We consider the orthonormal frame So = T, Ei, E 2 , E 3 
which is well defined everywhere in a neighborhood of the vertex p. Clearly the norm 
\U\, of an arbitrary tensorfield U, defined according to definition (27) coincides with 
the square root of the sum of squares of all the components of the tensor relative 
to this orthonormal frame. It is easy to se that, 

|A- (R*R)| < |A| • |R*R|. 



On the other hand, if e 4 = L, e 3 = L the null pair (68) and denote by a(R * 
R),/3(R * R), p(R * R), cr(R * R),/?(R * R), a(R * R) the null decomposition of 
R * R, as a Weyl field, relative to the null pair es, ca we can easily check that, 

|R*R| 2 % |a(R*R)| 2 + |/3(R*R)| 2 + |p(R*R)| 2 (88) 
+ |cr(R*R)| 2 + |/?(R*R)| 2 + |a(R*R)| 2 (89) 



13 Here $\ denotes the angular Laplace-Beltrami operator on the 2-surfaces St- 
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Indeed if we denote the Weyl field R * R by W and introduce its electric and 
magnetic parts Eij = Wiojo, Hij =* Wiojo we have, 



W\ 2 = \E\ 2 + \H[ 



Indeed, we have, 



Wijko — — &ij H s k, *Wijko —^ij E s k 

Wi 3 ki = - ^ijs^kit E st , *W l jki = - ^i 3S Gkit H st 

On the other hand, in terms of the null decomposition of W, relative to ei, e2, e 3 , e\, 

E a b = \&ab + \OLab ~ ^P&ab H ab = — | (l ab + \ a ab - ^a5 ab 

EaN = 2 §_a + 2@ a ^ aN = 2 * (La ~ 2 * @ a 



Hence, 



Enn — P Hnn — a 



\W\ 2 = \E\ 2 + \H\ 2 < \a\ 2 + |/3| 2 + \p\ 2 + \a\ 2 + \f3\ 2 + \a\ 



which proves (88). We now estimate the right hand side of (88). Clearly any 
null component of R * R can be expressed as a quadratic expression in the null 
components of R. We observe that no null component of R * R can be quadratic 
in a. This can be easily proved by a signature consideration. Indeed we assign 
signature 2 to a(R), signature 1 to /3(R) signature to p(R) and ct(R), signature 
— 1 to (3(R) and signature —2 to a(R). Similarly we assign signature 2 to a(R*R), 
signature 1 to /?(R*R) signature to p(R * R) and ct(R*R), signature —1 to 
(3(71*11) and signature —2 to a(R*R). It is easy to check that in the algebraic 
formula expresses the null components of R * R in terms of a quadratic form in 
the null components of R the total signature of each term must be the same as the 
signature of the corresponding null component of R * R. Thus, 

a(R*R) = (R*R) o3 63 = QrMp,a)]+Qr[M] 

where Qr[, ] denotes a simple quadratic expression in the null components of 1Z. 
Similarily, 

/3(R*R) = Qv[a,p]+Qv\l,(p,a)] 
p(K),a(R) = Qr[a, a] + Qr[/3, 0\ + Qr[(p, a), (p, a)] 
0(R*R) = Qr[/3, a] + Qr[(p, cr),/3]) 
a(R * R) = Qr[(p,a),a]+Qr\/3,0\) 

We now introduce the notation 

(\R\^) 2 = \a\ 2 + \p\ 2 + \p\ 2 + \a\ 2 + \P\ 2 (90) 

and deduce the following inequality, 

|R*R|<|R|t.|R| (91) 
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We are now ready to estimate the term I{p) . Using the bounds for lZ(p) of theorem 
4.13 and lZ(p, 5) < TZ(p) we derive, 



\m\ Z I IaiirMri 

JM-{p,S) 

< (/ I^A/ |A| 2 |R| 2 ) 1/2 

< TZ(p,8)( / ||R(i)||!oo||A(t)||| 2(St)( ii) 1/2 

Jt(p)-S 

< n{p,5)( ||A(t)||| 2(St) ) 1/2 sup ||R(t)|Ucc (St) 

Jt(v)-S t(v)-S<t<t(v) 



>t(p)-5 t(p)-5<t<t(p) 

We therefore have 

\I(P)\ < n(p,S)\\A\\ L2 ^- M) \\R\\ L ^- {piS)) (92) 

5.5. Estimate for J(p). In view of the fact that, 

/ \TL(;;L,L)\ 2 <K(p,S) 2 , 

JM- (p,«) 

we deduce, proceeding exctly as for /, 
Therefore, 

\J(P)\ < K( Pl 5) ||A|| L2(A ^- (Pi5)) \\R\\ L ~ W - M) (93) 

5.6. Estimates for L(p). We proceed as follows, 

\L(P)\ < I H|A||R|<(/ H) 1/2 (/ |A| 2 |R| 2 ) 1/2 

< llMlk 2 (7V-(p,5)) II A I|l 2 (^-(p,(5)) l|R|U~(A^-(p,5)) ( 94 ) 

5.7. The term K(p). Integrating by parts we rewrite K(p) as follows, 

K( P ) = - f fA-fK+f ( a T> a A ■ R 

JM-(v.6) JN-(v,S) 



We now estimate as follows, 



\K( P )\ < I (|)PA||y7R| + |C||y7A||R|) 

JN'-(v.S) 



IH- lp,5) 

< \\?A\\Lzw-( P ,5))(\\?R\\mAf-( P j)) + IICR||l 2 (^-( P ,5))) 

Therefore, 

\ K (P)\ £ liyA|| I ,2 (J v- (j)i5)) (||yR|| I ,2 ( ^- (P;5)) + ||R|U~(^-(p,5))IICIU 2 (Ar-( P ,5))) 
Going back to (86) and using the estimates for /, J, K, L obttained above we derive, 
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Proposition 5.8. The following estimate holds for all p G JA* and < S < 

|R-(p)l ^ £ + l|R-IU-(7v-(p,5))l|A|| z ,2 (JV -- (Pi5)) (7e(p,^) + ||/i||z,2 (A r-(p, 5 ))) 

+ \\fM\L2(Af-(p,8))(\\fR\\L2(Ar-(p,6)) + ||H-|U=»(7V-- IICII Z, 2 (AT- 

We now recall that according to (77) 

ll(p,6)<1l(p) <C(U,A )1Z . 

Furthermore, by (82) 

II^I|l 2 (A/'-(p,5)) < C(i*,Ao,7?o) 

and 



/ /-t(p) /■ 
\Jt(p)-6 J s t 



( ftip) \ 2 

< <5 sup / |C(*,^)| 2 

a- 2 \^*(p)-' 5 / 

< SC{U,A ,Tl ). 

Finally, in the next section we will establish that 

\\(t(p)-t)A\\ Lao{Ar - M) < C, \\fA\\ L 2 {// -(p,s)) < C 

with a constant C = C(i*, A ,1Z ). This in particular implies that 

l|A|| L 2 (Ar - (Pi5)) < 5?C(t*,A ,1l ). 

Putting this all together we deduce, 

l|R-(*)lk~ £ £ + 5^ sup HRIIioo^-^)) + sup ||yR||i,2 (Ar - (p>5)) 
pes* P eSt 

On the other hand, according to (80) we have, 

sup \\fK\\^(M-(p,5)) = sup.F«(p,<5)< ||DR(i-<S)|| L 2 + ( / ||R(t')|||=odt') 1/2 
Therefore, 



I|R(*)I|l~ < £ + \\im{t - 6)\\ L 2 + 6* sup ||R(i')lk- (95) 

t'e(t-s.t) 



5.9. Estimates for the error term £ . We first observe that, 

|nt|<i, |Dt|<i. 

Therefore, since /' and /" vanish for \t — t(p)\ < 5/2 and 

link- < <r\ ||/"|koo<r 2 
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we derive, 

\£\ £ ^||A|| L2(A ,- (P . 5) sup (r^lR^ll^csjjdt + ^HDR^IU^sj))^ 

t'e[t-s,ts/2] 

< sup (^IIR^IIl^^ + IIDR^IMs;)) 

t'e[ts,t-s/2] 

< sup (S- 1 (||R(t')IU 2 (E t ,) + ||DR(t')|U 2( E t ,)) + ||D 2 R(t')||L 2( E t ,)) 

t'e[t-5,t-5/2] 

The last step can be justified by a simple integration by parts argument. 

5.10. Final Estimate. Returning to (95), taking a supremum in t over an interval 
of size 8 and using the Sobolev inequality of corollary 2.10, 

||R(f)IU~ <C(\\R(t')\\ L i + ||DR(f + ||D 2 R(f)|U0 , 

we obtain, 

Proposition 5.11. There exists a positive S > 0, sufficiently small but depending 
only on A ,lZo,t*, such that the following estimate holds true, 

||R(t)IU~ < C6- 1 sup (||R(t')||L 2 + ||DR(t')||L 2 + ||D 2 R(t')||LO 

t'€[t-2S,t-5/2] 

with C a constant depending only on Ao, t*, IZq. 

We now return to propositions 3.7 and 3.8. Combining them with the proposition 
above we deduce, 

||DR(i)||£ a < C(||DR(t-5/2)||£ a + f HRCOIIioodf) 

Jt-S/2 

< OS" 1 sup (||R(t')||L 2 + ||DR(t')|li 2 + ||D 2 R(t')||| 2 ) 

t'€[t-6,t-5/2] 

ft 



||D 2 R(t)||| 2 < C(||D 2 R(t-«5)|| 2 2 + / ||DR(i')||| 2 ||R(t')ll!oo^) 

Jt-S/2 

< C\\T> 2 R(t-6)\\ 2 L2 



+ CS- 2 sup (||R(0l|L 2 + ||DR(t')Hi 2 + ||D 2 R(t')||| 2 ) 2 

t'e[t-5,t-S/2] 

Consequently, for some C depending only on /S.q,TZq and t*, 

Wnm^SCS- 1 sup ||R(t')ll^ (96) 
t'e[t-s,t-s/2] 

where, 

||R(i)||* a = ||R(t)||^ + l|DR(t)|| L2 + ||D 2 R(t)|| i2 (97) 
Iterating the estimate as many times as needed, in steps of size 6/2, we derive, 

Theorem 5.12. Assume that (M,g) is a globally hyperbolic extension of So ver- 
ifying the assumptions Al and A2. Let M* = U te [[_ t0;ti> )£ t C M with to = —1. 
There exists a constant C > depending only on A ,i* and initial data ||R(£o) I Iff 2 
such that, 

sup ||R(i)||if2 < C (98) 
te[t ,t.) 
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6. Proof of Main Theorem 1.1 

Theorem 5.12 established above provides us with global uniform bounds for the 
curvature tensor R and L 2 bounds for its first two covariant derivatives. Using 
elliptic estimates we can also derive L? bounds for the first three derivatives of the 
second fundamental form, see theorem 8.7. To finish the proof of the Main Theorem 
we only need to apply the following local existence result 

Proposition 6.1. Let (E*,<?, k) be initial data for the Einstein vacuum equations 
satisfying the constraint equations. We assume that E is compact and has constant 
mean curvature r = g^kij = const < 0. Let R denote the Ricci curvature tensor of 
g. Then there exists a smooth future Cauchy development of (E, <?, k) containing the 
region U te [ r T+p ]E t , where each E t is a constant mean curvature hypersurface (with 
mean curvature equal to t) and E T = E*. The constant p here depends only on the 
diameter and radius of injectivity o/E, the strictly negative constant r (t < 0) and 
the following constant, 

Tl* = PII^(E) + l|Vi?|| L2(s) + ||V 2 i?|| L2(s) (99) 

+ ||fc||L«(E) + l|Vfc|| L2(S ) + ||V 2 fc|| L2(s) + ||V 3 fc|| L2(s) . 

Proof : The proof requires a slight modification of the local exitence theorem 10.2.1 
in [C-K]. ■ 

Theorem 5.12 above, combined with the bounds on the second fundamental form 
stated in theorem 8.7, proved below in the Appendix, implies that for each hyper- 
surface E t C A4* with to < t < t* the constant defined in (99), is uniformly 
bounded. On the other hand, Theorem 2.9 together with L°° bounds on curvature 
implied by Theorem 5.12 guarantees a uniform bound for both the diameter and 
radius of injectivity of E t for to < t < i*. As a consequence, under assumptions 
Al, A2, as long as < we can construct a smooth globally hyperbolic CMC 
development containing the region U te [ t0iti> ]E t . 



7. Estimates for A. 

Proposition 7.1. Let A be the tensor defined in (84). Then for all < S < 

\\{t{p) -t)A|| L oc (Ar - (Pi5)) < C{U,A ,Tl ). (100) 

Proof . Recall, see definition 3.4, the convention b < 1 for an inequality of the 
form b < C(t*, A , TZo). We claim that it suffices to prove the proposittion for 
the case when A M is a vectorfield. The general case can be derived by a simple 
induction argument. Recall that we have, 

D L A + itr X A = 0. 
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with (sA) prescribed to be Jo at the vertex p. In view of the identity ^| = — (ntp) 1 
as well as estimate (67), it suffices to prove the inequality 

\\ s M\L°°(Af-(p,5)) ^ 1- 

Letting B = sA we have 

DiB = -l/2(H-trx)B, B| S=0 =J (101) 

Recall that, 

|B| 2 = B 2 + \B\ 2 = 2\B \ 2 + < B,B > 

where Bq =< B, T >, B_ is the projection of B on the foliation S t and < B, B >= 
g^B^B^. We shall first estimate < B,B > by observing that, 

^ < B,B >=-(?- tr X ) <B,B>, 



which in turn implies 
d 
d\ 

Therefore, since | < B,B > (0)| < |J(0)| < 1 



^ < B, B >= mp{- - trx) < B, B >, 
at s 



| <B(t),B(t) > | < | < B(0),B(0) > | exp ( / rup |tr X - -\dt') < 1, 

Jt s 

where the last inequality follows from (67) and (81). Therefore, for all t(p) — S < 
t < t(p), 

| <B(t),B(t) > | < 1. (102) 
We shall next derive a transport equation for B using the fact that B = — B a T + B_ 

= < D L B,T > + < B,D L T > 

as 



Observe that, 



= l/2(^-tr X )Bo+<S,D L T> 



< R, D L T >= -\<p- Y ( < B, D T T + DatT > ) 



Therefore, recalling our condition (28), 

I < B, D T T + DatT > | < A \B_\ 

Therefore, 

^|B | < |S | + A |S| 
On the other hand, from (102), | - B 2 + \B\ 2 \ = \ < B,B > | < 1 from which, 



\B\< ^/l + B 2 <\B Q \ + l 

Therefore, 

j t \B Q \< |B |+Ao(|S | + l) 
from which we deduce the estimate 

\B \ < 1 
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Thus, together with (102), we derive, 

|B| < |So| + \B\ < 1 (103) 
as desired. ■ 



7.2. Estimates for ^A. 

Proposition 7.3. Let A be the tensor defined in (84). Then, 

|| sup (t(p)-t)i\fA(t)\\\ Ll <l, (104) 

t(p)-5<t<t(p) 

||yA|| L 2 (Ar - (p , s)) < 1, (105) 

Proof In what follows we recall that points on Af~ (p, S) are parametrized by the 
coordinates (t,oj) with w£§ 2 . According to to Proposition 4.16 the volume forms 
dAtf-fag) on Af~(p, 5) and dSt on St are respectively equivalent to the expressions 
(t(p) — t) 2 dt da S 2 and (t(p) — t) 2 da S 2 with dcr§2 denoting the standard volume form 
on § 2 . Similarly equivalent are the quantities (t(p) — t) and s(t). Notation L 2 , above 
refers to the L 2 norm with respect to the measure dag2 . 

We begin by applying the results of lemma 4.9 to the equation DlA + ^tr%A = 0. 
and derive, 

FlWA) + Xabf b A^ = -^y a (tr X A M ) - i(C„ +2 o )trxA„ + R* Lo A A 
Therefore, 

fl(?a A ») + ^XifaK) = -Xabf b A^ - ^(f a tr X )^ - \ (Ca + ?7 a )tr X A„ + R„ A Lo A A 
which we rewrite in the form 

y7 L U aAI + tr X U aAI = -Xa h U fcAl + F QM , U| s=0 = (106) 

Fan = -\(fa^X)^ -\(Ca+ Vj^X^ + R^Lc^A (107) 

with U aAI = y a A M . Observe that 

|U| 2 = 2|U a0 | 2 + U aAl U^ 

Then, 

^ L U a0 + tr X U a0 = -Xa 6 U b0 + F a0 + ^^(n-^.n + fc;vj)U aj 
As a consequence, 

V7 L ( s 2 U a0 ) = .s 2 - tr X - x) ' U a0 + F a0 + ^^(n^n + k Nj )U a ^j := s 2 G 
Using that ^| = — n<yC and 

-^(s 4 |U a0 | 2 ) = s 4 g(G,U a0 ) 
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together with boundedness of n and tp, we estimate in the range t(p) — S < t < t(p), 

r t(p) 



sup S (t) 3 |U a0 | 2 |Ui < Hsup^t)- 1 / s(r) 4 g(G, V a0 )dT\\ 1 
t t Jt 

rt(j>) . i 
< Hsup^Ci)- 1 / S (x)S|G|^|| i2j HsupsCr^jUaol 2 !!^ 

t J t T 



< 



^Hsup^)- 1 / s(T)i\G\dr\\ 2 L2 +6||sup S (T) 3 |U a0 | 2 |Ui 

t Jt T 

To control || sup t s(t) _1 jf P ^ s(r) i \G\dr\\ L 2 we first estimate the integral, 



r {p> s 2 

J:= Hsup^i)- 1 / s(T^\(—-tr X -x)-Vao\dT\\ L 2 : 
t Jt S \ T ) 



as follows, 

r t(p) 



^< II sup / |(_-trx-x)|dr|| i »||sup*(r)'|U o0 ||| i j 



i /" t(p) 2 I 3 

< || sups(i)2 / |( _trx-x)| 2 dT||| ro ||sup S (T)^|U a0 ||U2 

t J t 8(T) r 

< (5^||sups(T)3|U a0 ||| L 2, 

r 

where the last inequality follows from (81) and (82). On the other hand, in view of 
(100), (81),(82) 

rt(p) 

Usurps^ 1 / s(T)^ftvx\\A\dT\\ Ll <S^\sn V (t(p)-T)\ftv X \\\Li<5K 

t Jt T 

MP) r PHP) 

Hsups(t)- 1 / 8(T)*\C + ri\\trx\\A\dT\\ Li < || sup(*(p) - r)|tr X | / |C + rfdrdu < 1. 

t Jt T Jt 

Moreover, using Proposition 4.13, we have 

s(r^\K^ La \\A\dr\\ Ll <\\ S ups(t)-^ jf (t(p) - r)|R M A io |dr|| L j 

^ 11^ La\\L 2 (Af-(p,8,)) ^^0- 

Note that it is the presence of an L component in the Riemann curvature tensor 
R^ A La which allows us to express it as a linear combination of the tangential terms 
a, (3, p, cr, (3 entering into the expression for the curvature flux. 

Finally, 

r t(p) 5 3 

Hsups(t)- 1 / aW'l^-^n - y.-n + fcATjOUajIdTllij < A <5|| sups(r)2 |U„.||| L a 

t Jt r 

Therefore, 

||sup S (i) 3 |U a0 | 2 || L i <e- 1 (l + S (t) 2 + ^o) + (e + e- 1 A 2 S (i) 2 )||sup S (t) 3 |U a .| 2 || Ll 

t r 

Combining this with the similar estimate on Ua^U^ we obtain 
|| sup (t(p)-tf\U a0 \ 2 \\ Ll <1 + K , 

*(p)-<5<t<*(p) 
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which gives (104). The argument above also provides the inequality 

1 f t{p) 

(t(p)-t)\\U\\ Li < ( t (p)-t) 3 + -^_ jf ( (t (p) -r)||R^ L JUj + ||C + !?lk) 
Using the maximal function estimate we then obtain 

l|U|| L2(A A- (pA ))<l + ftJ 
and hence (105). ■ 



8. Appendix 



Recall that the curvature tensor R can be decomposed into its electric and magnetic 
parts E, H as follows, 

E(X, Y) =< R(X, T)T, Y >, H(X, Y) =< *H(X, T)T, F > 

(108) 

with *R the Hodge dual of R. One can easily check that E and H are tangent, 
traceless 2-tensors, to S t and that |R| 2 = \E\ 2 + \H\ 2 . We easily check the formulas 
relative to an orthonormal frame eo = T, ei, e2, e^, 

RafecO = — Gabs H sc , *R a fccO =^abs E sc (109) 

Rafecd = ^abs^cdt E s t, *Ra6cd = — ^abs^cdt Hst 

We recall below some of the main formulas involving k, E and H. 

Eij-Rij = tvkkij -kfk 8j (110) 
Hij = curl kij (HI) 
where, for any given symmetric two tensor I of S t one defines 

curl hj ^ef VJ bj + ef VJ ib . 

We also recall the constraint equation for k, 

V j kij - Vitrfc = (112) 

In the particular case when trfc is constant equations (111) and (112) form an elliptic 
Hodge system on E t , 

div k = curl k = H (113) 

8.1. Elliptic L 2 - estimates for Hodge systems. Here we recall the following 
lemma concerning rank-2 symmetric Hodge systems on a 3 dimensional compact 
Riemannian manifold S. 

Lemma 8.2. The following L 2 elliptic estimates hold on a 3 dimensional Riemann- 
ian manoflod S. 



i. Let V be a symmetric tracelss 2— tensor on S verifying, 

div V = p, curl V = a 



(114) 
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Then, 

£ (|W| 2 + m mn V lm V t n \r\V\ 2 ) = Jj\<j\ 2 + ±\p\ 2 ) (115) 

where Rij is the Ricci curvature of T, and R its scalar curvature. 
ii. For a scalar <f> we have, 

[ |V 2 ^| 2 + f WVrfVrf = f |A0| 2 
Proof . See Proposition 4.4.1 in [C-K]. ■ 



8.3. Apriori estimates for k. We now apply lemma 8.2 to the Hodge system 
(113) for k on a fixed hypersurface S = S t , t < 0, 

£ (|Vfc| 2 + 3(k 2 ) mn (E mn + (k 2 ) mn ) - i|fc| 2 |fc| 2 ) = £ \H\ 2 

Interpreting k as 3 x 3 symetric matrices we can write 

i{k 2 ) mn {k 2 ) mn - \\k\ 2 \k\ 2 = 3tr(fc 4 ) - i(trfc 2 ) 2 

Observe that we have the pointwise inequality , for an arbitrary symmetric matrix 
k, tr(fc 4 ) > ±|fc| 4 . Therefore, 

I |Vfc| 2 + i|fc| 4 < jf \H\ 2 + £ \E\ \k\ 2 < jf |ff| 2 + |£| 2 + ^ J|fc| 4 
This proves the following: 

Proposition 8.4. On any leaf £ of a constant mean curvature foliation S t the 
second fundamental form k verifies the estimate, 

jf |Vfc| 2 + J|fc| 4 < ^ |ff| 2 + |£| 2 = |R| 2 . (116) 
In view of the energy estimate of proposition 3.3 we derive, 

Corollary 8.5. The following estimates hold true with a constant C depending 
only on Ao and t*, 

\\Vk(t)\\ L 2 + \\k(t)\\ L ,<CTZ (117) 



Indeed diagonalizing k it suffices to prove the inequality for arbitrary real numbers a, b, c, 
3( a 2 +5 2 +c 2 )2 > a 4 +5 4 +c 4^ 



CURVATURE 



37 



8.6. Higher derivatives estimates for k. To derive second derivative estimates 
for k we rewrite the curl equation in (113) in the form, 

Differentiating we obtain, 

or, symbolically, 

Ak = R*k + VH 

where R * k is a quadratic expression with respect to the Ricci curvature R of S t 
and k. Thus, since, the Ricci curvature R can be expressed in the form, 

Rij ki a k j -\- tikkij — Eij 

we derive, 

|Afc| < |fc| 3 + |^||fc| + |Vi/| 

Therefore, 

/ \Ak\ 2 < f (|fc| 6 + |£| 2 |fc| 2 + |Vtf| 2 ) 
It is easy to see by a standard integration by parts argument that, 
/ \V 2 k(t)\ 2 < f \Ak\ 2 + f \R\ 2 \Vk\ 2 

Consequently, 

j |V 2 fc| 2 < [ (\k\ 6 + \E\ 2 \k\ 2 + \VH\ 2 ) + [ \R\ 2 \Vk\ 2 

Therefore, since < A , 

\\V 2 k{t)\\ L , < A 2 (||fc(t)||l 4 + ||R|| 2 2 ) + ||V^(i)|t 2 2 (118) 

+ \\R(t)\\ L ~\\Vk(t)\\ 2 L i 

It is easy to see that \\VH(t)\\ 2 L2 < ||DR(t)|| L 2. Also, 

\\R(t)\\L- < IRII^ + llfclll^^HRWII^+A 2 

Therefore, in view of theorem 5.12 and the bounds for ||Vfc(£)||i2 and ||A;(t)||£4 of 
corollary 8.5 we derive from (119), 

||V 2 fc(t)|| L2 < C 

with C a constant depending only on A ,i* and TZq. 

Differentiating once more the equation for Afc and proceeding in the same fashion 
we can also derive similar bounds for the third derivatives of k. This proves the 
following. 

Theorem 8.7. The second fundamental form k of the t foliation satisfies the fol- 
lowing estimate, for all to < t < t*, 

||V 3 fc(t)|| L2 + ||V 2 fc(t)|| i2 + ||Vfc(t)|| L2 + \\k{t)\\ L < < C (119) 

with C a constant depending only on A ,i* and TZq- 
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